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INTROuUCTlUN 

Part  1  of  this  report  indicated  a  novel  solution  to  tne  basic  diffusion 
equation  of  Pnysics  where  the  field  ooundary  extends  from  zero  to  positive 
infinity.  The  nodal  points  of  tne  field  net  are  identified  as  terminating 
polynomials  with  the  numerators  of  the  coefficients  found  first  by  deduction  - 
for  the  lower  orders  -  and  then  by  extrapolation.  Part  11  considers  the 
numerical  analysis  employed  to  complete  tne  entire  set  of  tables. 

PROCEDURE 


A.  Direct  Differences 

Formulation  of  the  polynomial  form  of  the  discrete  solutions,  Eq.  (1), 
of  the  diffusion  equation  from  the  Schmidt  plot  geometry  is  described  in 
References  (1)  and  (2)  and  reflects  a  progressive  trigonometric  construction 
where  the  degree  and  term  extension  increases  with  time  and  decreases  with 
distance,  time  and  distance  referring  to  the  unsteady  heat  flow  application. 


T  (N.P) 


-r  [  A  -  Aj^m  +  A  -  A^  m^  +  A^m^ -/••••+ 


where  m  is  the  independent  variable, 


N  is  a  distance  index, 

P  is  a  time  index, 

_  (Pi-iO-2-  I  sin  (P+iOT^/Z  I 

2  ■  . 

j  =  (k,  -  term  exponent  of  m) ,  the  individual  term  denominator  exponent, 
and  y  =  m  (T^  -  T^U),  witn 

A,  A,,  . Aj^  the  numerical  coefficients  of  the  interior 

terms  of  the  equation. 

The  numerators  of  each  term  of  the  "nodal"  equations,  are  uniquely  related 
to  adjacent  time  and  distance  term  coefficients  of  the  same  degrees.  This 
relation,  originally  found  accidentally,  is  correlated  by  the  table  of 
differences  shown  in  Table  1  &  1-A  where  the  boxed  vertical  sequence;  17548, 
25147,  35401,  49024  and  66868,  is  established  by  the  reduction  of  the  Schmidt 
plot  through  the  trigonometric  analysis.  Step-wise  right  moving  subtraction 
generates  a  column  of  residual  zeroes,  an  adjoining  column  of  ones,  and  a 
digital  sequence  identified  as  "IV"  in  Table  1.  Reducing  this  column  "IV"  to 
zero  vertically  then  allows  a  corresponding  determination  of  the  particular 
values  of  the  entire  matrix  from  inspection  of  the  biased  rows  I,  II  and  III. 
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Table  1-a  (continuation  of  Table  1) 
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B.  Suniaing  Progression 


Rewriting  Table  1  as  Table  2,  where  the  biased  rows  are  horizontal  and 
the  zero  colunn  is  left-justified  instead  of  right-justified,  shows  that  the 
suu  of  any  two  adjacent  column  values  of  any  row  gives  the  value  of  the  next 
row  entry  directly  under  the  right-wise  addendum.  The  sequence  of  Row  1, 
regardless  of  the  degree  of  the  term  represented  by  the  matrix,  always  starts 
with  zero  and  then  maintains  alternate  zeroes  to  infinity.  Each  of  the 
difference  tables  corresponding  to  a  given  "m"  exponent  (Eq.  (1))  can  be 
resolved  similarly,  except  that  each  first  row  is  carried  in  a  unique 
progression.  Tnis  progression  is  repeated  without  the  Interspersed  zeroes 
within  the  matrix. 


C.  Pascal  Triangles 

The  classical  Pascal  triangle.  Table  3,  can  be  formed  oy  simple 
addition  where  each  term  is  the  sum  of  the  two  previous  superior 
terms  and  individual  entries  are  represented  by  the  binomial  coefficient, 

2  Z  ^ 

(  )  =  where  z  and  w  represent  the  row  and  colum  of  a  particular 

coeYf  icient  of'^a'  binomial  expansion.  A  modification  of  the  Pascal  triangle  is 
found  by  writing  the  diagonals  as  rows  which  then  generates  the  "arithmetic 
square",  Table  4,  also  known  historically. 


It  is  precisely  these  progressions  alternating  with  zeroes,  which 
comprise  the  first  rows  of  the  individual  "summing  progressions"  shown  as 
Table  2.  In  this  case  the  digital  enumeration  of  the  rows  indicates  the 
degree  of  the  "m"  tern. 

D.  Pockhammer's  Symbol 


Each  row  of  Table  4  can  be  examined  by  finite  differencing  to  establish, 
via  Gregory-Newton,  a  definitive  polynomial  expansion  extending  to  infinity. 
Furthermore,  the  algebraic  equation  can  be  simplified  to  a  factorial  form 
known  as  Pockhammer's  Symbol  or  as  a  n  factorial.  Appendix  A  presents  an 
example  of  such  a  uevelopment  for  "m"  degrees  zero  through  three.  The 
complete  arithmetic  square.  Table  2,  can  then  be  written  as 


f(v) 


where 


1 

(r+D! 


(2) 


r  =  degree  of  m, 

V  =  column  value 

and  f(v)  “  row  value  of  Table  2.  From  Table  2,  a  complete  construction  of 
Table  1  follows. 


Table  2.  Summing  Progression  for  Degree  zero  Numerators 


Table  3.  Pascal  Triangle 
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Sunimary 


A  straight-forward  aiethod  (aritnmetlc  squares)  is  described  to  permit  the 
numerical  construction  of  the  differencing  tables  of  Part  I  of  this  report. 

The  derivation  through  the  Pascal  triangle  and  correspondence  to  Pockhammer's 
Symbol  notation  is  demonstrated. 


lit 
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APPENDIX  A 


Determination  of  Pockhammer's  Notation 

A  discussion  of  the  Gregory-Ne»7ton  analysis  is  presented  in  Reference  4.  It 
consists  of  determining  a  polynomial  expression  to  represent  a  progressive 
sequence  of  numbers.  In  tiie  present  application,  it  is  used  to  examine  the 
base  row  development  of  Table  4. 

Given  a  unit  stepping  difference  in  a  counting  reference,  v,  and  a  matched 
sequence  f(y); 


f(y)= 


f(v)  + 


Af(v) 

1! 

3, 


„(1)  ^  A^f(v) 
^  2! 


^(v)  (3)^ 

3 !  ^ 


(2) 

y 


wnere 

V  is  tne  step  level, 

f(y)  is  the  dependent  variable, 

4^f(v)  are  Liie  diagonal  values  of  the  difference  table,  and 

y^^^  =  1 

y(l)  =  y 
y(2)  ^  y(y-l) 

y(3)  ,,  .y_2)  •••• 

etc . 

For  the  case  of  "m"  degree  aero  where  f(y)  is  from  Table  4: 
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f(y)  f(v) 


0 

1 

1 

2 

2 

3 

3 

4 

4 

5 

5 

6 

6 

7 

7 

8 

8 

9 

9 

10 

10 

11 

11 

12 

Af  (v) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


f(y)  =  f(v)  + 


Af(v)  ,(1), 


1! 


=  1  +  V  +  0 

=  (v  +  1) 
For  m  degree  1  : 

V  f  (y) 

0  1 

1  3 

2  6 


3  10 

4  15 


f(v) 

2 

3 

4 

5 


Af  (v) 


A^f (v) 


0 

0 

0 


f(y) 


f(v:  + 


Af(v)  ^(1) 
1!  ^ 


A^f (v) 
2! 
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Af(v) 


«  1  +  2v  +  -I-  v(v-l)  +  0  + 
*  -j  (v^  +  3v  +  2) 


*  Y  (v  +  1)  (v  +  2) 
For  m  degree  2 
V  f(y)  f(v) 


3 

6 

10 

15 

21 


3 

4 

5 

6 


A^f(v) 

1 

1 

1 


=  1  +  3v  +  y  (v)(v-l)  +  (v)  (v-1)  (v-2)  +  0  • 

=  1  +  3v  +  "I  (v^  “  ^  -  3v  +  2) 

=  4"  (v^  +  6v^  +  llv  +  6) 

6 

=  (v  -t-  1)  (v  +  2)  (v  +  3) 

The  pattern  continues  so  that  : 

Degree  of  "in" 

=  r  f(y) 

0  (v  +  1) 

1  J_  (v  +  1)  (v  +  2) 

2 

2  J_  (v  +  1)  (v  +  2)  (v  +  3) 

6 

3  _1_  (v  +  1)  (v  +  2)  (v  +  3)  (v  +  4) 

24 

4  1  (v  +  1)  (v  +  2)  (v  +  3)  (v  +  4)  (v  +  5) 
120 


A^f(v) 

0 

0 
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5 


1  (v  +  1)  (v  +  2)  (v  +  3)  (v  +  4)  (v  +  5)  (v  +  6) 
720 


and  ttie  general  expression  is 


1 

f(y)  - -  (v) 

(r  +  1)  !  r  +  1 

where  (v)  =  (v  +  1)  (v  +  2)  (v  +  3)  .  .  .  .  •  (v  +  r  +  1)  . 

r  +  1 


With  respect  to  the  original  time  index,  P,  of  the  diffusion  equation  polynomial 
V  =  P-1  by  Table  3,  and 

(v)  =  (P  -  1  +  1)  (P  -  1  +  2)  .  .  .  .  (P  +  r) 

r  +  1 

=  P  (P  +  1)  (P  +  2)  (P  +  3)  .  .  .  .(P  +  r) 

so  that 

1 

f  (y)  =  -  (p) 

(r  +  1)  :  r 

wnich  is  known  as  Pockha.umer 's  Syraaol. 


if 

Gravio  A.  Korn,  iiathematical  riandbook  for  Scientists  and  Engineers »  McGraw- 
Hill  Book  Co.,  Inc.,  New  York,  etc.,  1961. 
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List  of  Symbols 


h 

j 

k. 

m 

r 

V 

f(v) 

w 


z 

^  ‘  A3 ,  .  .  .  . 


U 

P 

T 

9 


exponent  of  2  in  external  denominator 
exponent  of  2  in  each  term  denominator 
exponent  of  "m"  in  final  term 
independent  variable 
degree  of  "m” 

column  value  of  stepping  sequence 

row  value  of  stepping  sequence 

inferior  component  of  binomial  coefficient 

superior  component  of  binomial  coefficient 

numerical  coefficients 

distance  index 

time  index 

dependent  variable 

external  numerator 
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